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The origin of the micro-Gauss magnetic fields in galaxy clusters is one of the outstanding prob¬ 
lem of modern cosmology. We have performed three-dimensional particle-in-cell simulations of the 
nonrelativistic Weibel instability in an electron-proton plasma, in conditions typical of cosmological 
shocks. These simulations indicate that cluster fields could have been produced by shocks propa¬ 
gating through the intergalactic medium during the formation of large-scale structure or by shocks 
within the cluster. The strengths of the shock-generated fields range from tens of nano-Gauss in 
the intercluster medium to a few micro-Gauss inside galaxy clusters. 
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The origin of the micro-Gauss magnetic fields observed 
in galaxy clusters ULUS poses one of the most intriguing 
problems in modern cosmology. The most common ex¬ 
planation invokes the amplification of seed or primordial 
fields by hydrodynamic turbulence that have been excited 
during the processes of large-scale-structure (LSS) for¬ 
mation. Although there are viable astrophysical mecha¬ 
nisms that can generate seed fields with B ~ 10" 16 Gauss 
or weaker [jj, 30- 0, recent cosmological simulations 0 
show that structure formation can amplify the field by 
no more than a factor of 10 3 . Hence, in order to explain 
the observed intergalactic field, one needs seed fields as 
strong as B ~ 10~ 9 Gauss. Scenarios with galactic winds 
and quasar-driven outflows mil can provide fields of 
such strength, but they are rather localized. It is thus 
not clear whether they can explain entirely the origin of 
the intergalactic fields in galaxy clusters. 

Here we show that magnetic fields can be produced by 
collisionless shocks in galaxy clusters and in the interclus¬ 
ter medium (ICM) during LSS formation. Cosmological 
iV-body and hydrodynamic simulations of LSS formation 
11,113 have shown that shocks with Mach numbers up 
to M ~ 100 are ubiquitous on scales of few to few tens 
of megaparsecs. Theoretical analysis of non-magnetized 
collisionless shocks indicates that they can generate sub- 
equipartition fields 00 , 0 . We verify this prediction 
with state-of-the-art numerical simulations. We present 
here three-dimensional (3D) particle-in-cell (PIC) simu¬ 
lations of the nonrelativistic (with v = 0.1c) Weibel insta¬ 
bility 00 in an electron-proton (with m p /m e = 100) 
plasma, thus guaranteeing a clear separation of the rele¬ 
vant time scales. These simulations are computationally 
expensive and represent a significant advance over previ¬ 
ous studies, which simulated relativistic shocks (v ~ c) in 
an electron-positron or low-mass-ratio electron-ion plas¬ 
mas (m,/m e < 16) 000. Note that a recently 
discussed possibility that cluster shocks may produce the 
magnetic fields seen in galaxy clusters 0 was based on 
the assumption that the results of relativistic simulations 
will also apply in the nonrelativistic regime. Our nonrel¬ 


ativistic simulations fully confirm theoretical predictions 
and indicate that LSS shocks can produce magnetic fields 
of strengths of tens of nano-Gauss to few micro-Gauss in 
the intergalactic medium (IGM) and ICM, respectively. 

The mechanism of the field generation at shocks is 
rather simple 0 0 - As a shock propagates into an 
ambient medium, it reflects (or scatters) a fraction of the 
incoming (in the shock frame) particles back into the 
upstream region which then form counter-propagating 
streams. Both groups of particles (ICM/IGM and re¬ 
flected particles) have bulk velocities of order the shock 
velocity w s h; they can also have some thermal spread. 
Both protons and electrons form the streams, so both 
species participate in the instability [32I . One can con¬ 
sider each charged particle in these streams as an ele¬ 
mentary current. Since like currents attract each other, 
it is energetically favorable for the elementary currents to 
merge into larger current filaments. This process is inhib¬ 
ited at scales smaller than the plasma skin depth, ~ c/oj p 
(w p is the plasma frequency), by strong electrostatic re¬ 
pulsion of like charges. In contrast, on large scales, the 
currents are quasi-neutral because of Debye shielding in 
a plasma. Hence, the filaments and associated magnetic 
fields grow rapidly. The process stops when most of the 
particles become trapped in the produced fields and can 
no longer amplify the field. This happens when the par¬ 
ticle Larmor radius pl = v_\_b/vj c (vj _b is the particle ve¬ 
locity component transverse to the local magnetic field, 
and u> c = eB/mc is the cyclotron frequency) becomes 
comparable to (or less than) the characteristic correla¬ 
tion scale \b of the field, pl/^b ~ 1. At this time, 
the particle distribution is effectively isotropized, and so 

^thermal ~ ^_L B ~ ^sh- 

The anisotropy of the particle distribution near a shock 
can be parameterized as, 

A = (e„ - ej_)/'etot =* (M 2 - 1 )/(M 2 + 1), (1) 

where cm oc V sh is the energy of the particle along the 
shock propagation direction; e± oc r’t hermal — c s is the 
thermal energy in the plane of the shock; e to t = en + ej_ 
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is the total energy; c s is the sound speed upstream; and 
the Mach number of the shock is M — v s h /c s - For strong 
shocks, M 1, the anisotropy parameter is close to 
unity, A ~ 1 . At a shock, the bulk velocities of the elec¬ 
tron and proton components are both comparable to the 
shock velocity. Hence, the protons dominate over the 
electrons in the overall energy budget, and the magnetic 
field generated by the electrons is negligible compared 
with that generated by the protons. The growth rate 
and the wavenumber of the fastest growing mode (which, 
in fact, sets the spatial correlation scale of the produced 
field) are 7 B = Auj PtP (v s h/c) and k B = Auj PtP /c, where 

u p,p = ( 47 re 2 n p /m p ) 1 ^ 2 « 1.32 x 10 z nJ 2 s -1 is the pro¬ 
ton plasma frequency, and n p and m p are the number 
density and the mass of the protons, respectively. (We 
use cgs units throughout, unless stated otherwise.) The 
kinetic calculation of the growth rate and the instability 
threshold have been examined elsewhere [22j. Order-of- 
magnitude estimates of the magnetic-field e-folding time 
and the held correlation length at strong shocks (M 1 ) 
are readily obtained as 

t b ~ 1 / 7 b - 2 x 10 2 v-^ >7 rijcM ,_4 s > ( 2 ) 

\ B ~ 2tt /k B ~ 10 10 njcM,-4 cm > ( 3 ) 

for a typical ICM proton density of mcM ~ 10 ~ 4 cm -3 
and a typical shock velocity v s h ~ 10 ' cm s _1 ; as 
usual, we denote nicM ,~4 = ^icm/( 10 -4 cm -3 ) and 
^sh ,7 = W(10 7 cm s _1 ). Since it takes N ~ few x 10 
e-foldings to produce strong fields, we can readily es¬ 
timate the thickness of a region of the held growth as 
A ~ Nt b i> sh ~ N X B . 

The saturation level of the magnetic held is estimated 
from As ~ Pl = 'Csh/w CjP , where u C:P = eB/rripC « 
9.58 x 10 3 B s ” 1 is the proton cyclotron frequency. In 
a multiple-species plasma, however, saturation occurs at 
equipartition with the lightest species 0]. To incorpo¬ 
rate this, we introduce an efficiency factor 77 , which in 
the electron-proton plasma is of order m e /m p . Finally, 


B 2 / 8tt 

m P n P V sh / 2 


B 2 

87rp sh 


A 2 p ~ 10 “ 3 , 


(4) 


where p s h is the gas pressure behind the shock, and the 
last estimate is for strong shocks, A ~ 1. 

Although there is no doubt that magnetic helds are 
generated at shocks through the Weibel instability, it 
is not clear whether they survive sufficiently far down¬ 
stream to produce longstanding magnetic fields. The 
concern arises from the fact that the wavelength of 
the fastest-growing mode in the linear Weibel-instability 
analysis is very small, X B ~ 2nc/u> Pt p ~ 10 10 cm for 
a typical ICM particle density of n ~ 10 -4 cm -3 . 
Therefore, it is possible that the extremely short spa¬ 
tial scales- i.e., sharp field gradients—can be rapidly de¬ 
stroyed by dissipation on a plasma time scale of t b ~ 


10 2 s. Should this happen, the fields would occupy only a 
very narrow region near the shock front and, thus, would 
not result in long-lived cluster helds. Recently, it has 
been shown 0| 1 both theoretically and using relativistic 
PIC simulations that the correlation length of the mag¬ 
netic held (and, hence, its gradient scale) grows rapidly 
with time, thus drastically reducing diffusive (Ohmic) 
dissipation. These results suggest that the magnetic 
helds produced should survive on cosmological times. 

In general, it is far from clear that nonrelativistic 
shocks can generate helds in the way relativistic shocks 
do. In order to test this, we have performed a set of 3D 
and 2D PIC simulations [13,0] using a state-of-the-art, 
massively-parallel, electromagnetic, fully-relativistic, 3D 
PIC code OSIRIS 2.0 [ijJ. In our PIC simulations, the 
initial conditions are taken to be two streams of electrons 
and ions moving with relative bulk velocity v s h , which in 
our simulations we take to be 0.1c. The four “species” 
of particles (the upstream and downstream electrons and 
ions) is then each assigned a Maxwell-Boltzmann distri¬ 
bution of velocities about the bulk velocity. Our ions are 
“light protons,” positively-charged particles with mass 
mions = 100 m e , a mass ratio large enough to guarantee 
that the electron and ion time scales are clearly sepa¬ 
rated. All of our simulations have volumes of 128 3 cells, 
although the cell sizes differ. 

We hrst ran four shorter-duration test simulations, 
three of them 3D and one 2D. In the first 3D run, the box 
size was (25.6 c/w Pje ) 3 (where w Pj6 is the electron plasma 
frequency), and there were 4 particles/species/cell. The 
second simulation differed from the first only in that 
the box size was (12.8 c/u; Pi e) 3 , and the third simula¬ 
tion differed from the hrst only in that it had 8 par¬ 
ticles/species/cell. These test simulations showed that 
neither the box size, nor the number of particles per cell 
affect the results within the available computational re¬ 
sources. We also ran a 2D simulation with 1280 2 cells, i.e, 
with a box of size (256.0 c/w Pje ) 2 , with periodic boundary 
conditions and 9 particles/cell/species in order to exam¬ 
ine the dynamics of a similar system in the plane trans¬ 
verse to the bulk motion of the shocked plasma (the con¬ 
figuration is identical to that in simulations of Ref. 0[). 
Comparison between the 3D simulations and the 2D sim¬ 
ulations do not show significant differences (e.g., in terms 
of the saturation level of the magnetic held), only re¬ 
vealing the limitations of the 2D simulation (the Weibel 
instability is stronger in 2D), and confirming that the 
transverse dimensions of the 3D simulation box are not 
strongly affecting the held dynamics on the time scales 
analyzed here. 

We then ran very long 3D simulations of colliding 
plasma slabs, for four sets of plasma parameters. One 
of the plasma slabs describes a shocked high Mach- 
number plasma (M = 20, Uth.eshock /c = 0.05, with ei¬ 
ther u t h,i shock/c = 0.005 or i>th,i shock/c = 0.0) with bulk 
motion along the X\ direction. The different ion ther- 
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mal velocities correspond to the two extreme cases of a 
strongly turbulent shock |28| , where electrons and ions in 
the shocked plasma are thermalized by plasma turbulence 
and a laminar shock {2lJ [30| ■ The second plasma slab de¬ 
scribes the IGM/ICM plasma with either cold electrons, 
T e = 0, or hot electrons, T e = Ti ~ 10 eV, such that 
vth,eiGM/c = 0.05, uth,iiGM/c = 0.005. All four of these 
simulations used a box of size (25.6 c/co P}e ) 3 . We found 
that the different physical parameters did not reveal any 
significant differences in the evolution of ts- Our present 
choice of simulation parameters is strongly limited by the 
time scales involved in the mechanisms described here, 
and it aims to illustrate the key features of the magnetic- 
field generation in conditions relevant for nonrelativistic 
collisionless shocks. 

In fact, when examining the temporal evolution of 
es measured in the 3D simulations (Figure QJ, we ob¬ 
serve the key role played by the ions, with most of the 
magnetic-field energy generated by the Weibel instabil¬ 
ity originating in the shocked ions. All the runs re¬ 
vealed cb — 10” 3 . Note that the Weibel-field growth 
of the ions saturates at lower relative eg than for elec¬ 
trons eB,i ~ (m e /mj) 1 / 2 es ie , where for species s = i,e, 
cb,s = (B 2 /8n)(m s n s v 2 h /2)~ 1 . A strong thermalization 
between the electrons in the two slabs is achieved very 
early in time via the electron Weibel instability, but ion 
thermalization is not observed in our simulations. Other 
instabilities with longer time scales (e.g., the ion acoustic 
instability) will be responsible for this. These instabilities 
are not observed in our simulations since the simulation 
box is not large enough. 

The time scale for energy transfer between the ions 
and the magnetic field is the time scale for the Weibel 
instability of the ions, longer than the electron Weibel 
instability by a factor of {m.i/m e ) 1 ^ 2 , thus making its ob¬ 
servation in numerical simulations very time consuming. 
The structure of the generated magnetic field depicted 
in Figure [3 shows the typical configuration of a Weibel- 
driven field, in 3D and in the 2D plane transverse to the 
bulk motion of the shocked plasma, surrounding the self¬ 
generated current filaments, which are already evolving 
to longer wavelengths. 

We have demonstrated that magnetic fields are pro¬ 
duced at nonrelativistic collisionless shocks and their 
strengths are comparable to that observed in clusters. 
It is thus natural to explain the observed fields by the 
Weibel instability. If so, then the magnetization of clus¬ 
ters should begin around the reionization epoch, at red- 
shifts of z ~ 10 — 20, when the gas becomes highly ionized 
and particle collisions become rare and inefficient. Our 
studies reveal that the magnetic field grows to an en¬ 
ergy density of roughly a tenth of a percent of the initial 
kinetic-energy density, and hence constitutes a similar 
fraction, es ~ 10~ 3 , of the thermal energy density of 
the shocked gas. The actual number depends on com¬ 
plicated nonlinear dynamics of the currents in the down¬ 


stream region. This value of the equipartition parameter 
corresponds to a magnetic-field strength of order 

B ~ 1(T 8 e ^ /2 _3 v sh ,7 ^icm ,—4 Gauss. (5) 

These values correspond to tens of nano-Gauss in the 
ICM and a few micro-Gauss inside galaxy clusters, the 
latter in excellent agreement with observations. 

The simulations presented here model a strong shock 
with Mach number, M = 20. Can a weak, M ~ 1, shock 
generate fields as well? 3D simulations of weak shocks are 
hardly possible at present. However, theoretical analysis 
of the Weibel instability shows that fields are generated 
when the shock velocity is larger than the thermal ve¬ 
locity of ICM/IGM particles by a factor of two or more; 
that is, when M > 2. The exact number depends on the 
actual particle distribution at the shock. 

LSS shocks can be observed via (i) synchrotron emis¬ 
sion by the shock-accelerated electrons in the in situ gen¬ 
erated magnetic fields; (ii) inverse-Compton scattering 
of cosmic microwave background photons by the shock- 
accelerated electrons; (iii) Sunyaev-Zeldovich effect on 
thermal electrons in the shocked medium downstream; 
and/or (iv) an abrupt change in the Faraday rotation 
measure across the shock. The shock front appears to be 
very thin and will likely be unresolved. Since no sign of 
proton thermalization is seen by the end of the simula¬ 
tions, t ~ 500^^, we can put a constraint on the shock 
thickness, 

A sh > 10 13 /j 2 _ 4 cm. (6) 

Our present analysis does not consider the evolution 
of the fields on cosmological time scales. Numerical sim¬ 
ulations of this type are hardly possible within the next 
few years. Theoretical considerations suggest that in¬ 
verse cascade should result in the rapid transfer of mag¬ 
netic energy from small (shock) to large (cosmological) 
scales, thus leading to the long-term survival of the fields 
241. An alternative possibility, that other instabilities 
M jresent in plasma could destroy the field entirely, is 
unlikely because the time scales involved are long enough 
for dynamos driven by turbulence and sheared motions 
of gas in clusters to further amplify and preserve the 
shock-generated magnetic fields. A detailed study of 
these issues is highly important, yet extremely difficult 
and should be addressed in the future. 
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FIG. 1: The evolution of the magnetic field energy normalized 
by the total initial kinetic energy, eg, is shown with the light 
blue line. The energy in the magnetic field is predominantly 
associated with the field components parallel to the shock 
plane. For comparison, the similarly normalized energies for 
the four particle species are also shown. 



FIG. 2: Magnetic field energy density at t = 2000/uJ p , e (u) P ,e 
is the electron plasma frequency). The blue iso-surfaces corre¬ 
spond to a value of eg ~ 8x 1CU 3 . The projection in the X2-X3 
plane (the shock plane) is the value of eg averaged along xi 
(the shock propagation direction) with red color correspond¬ 
ing to a peak value of eg ~ 6 x 10 -2 . The color scale in the 
projection plane is linear. 


do not consider the Langmuir instability in this paper. 










